Charge transport in purple membrane monolayers: A sequential tunneling approach 
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Current voltage (1-V) characteristics in proteins can be sensitive to conformational change induced 
by an external stimulus (photon, odour, etc.). This sensitivity can be used in medical and industrial 
applications besides shedding new hght in the microscopic structure of biological materials. Here, 
we show that a sequential tunnehng model of carrier transfer between neighbouring amino-acids in a 
single protein can be the basic mechanism responsible of the electrical properties measured in a wide 
range of applied potentials. We also show that such a strict correlation between the protein structure 
and the electrical response can lead to a new generation of nanobiosensors that mimic the sensorial 
activity of living species. To demonstrate the potential usefulness of protein electrical properties, 
we provide a microscopic interpretation of recent I-V experiments carried out in bacteriorhodopsin 
at a nanoscale length. 

PACS numbers: 87.14.E-; 87.15.hp; 82.39.Jn; 33.15.Hp 
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The nature of electron transport (ET) in biological ma- 
terial is matter of outstanding interest both for the pure 
speculative aspects and for the applications. As a mat- 
ter of fact, organic/biological-based devices are the new 
frontier of technology, due to their potential low cost, low 
size, high specificity, etc. [H, Q] . Usually biological mate- 
rials are not easy to be investigated, because a standard 
way of preparation is still not available Neverthe- 
less, there are some relevant exceptions, like monolayers 
of purple membrane (PM), a part of the cell membrane 
of the halophile Halobacterium salinarum, which is easy 
to prepare and suitable for direct measurements. PM is 
constituted by a single type of protein, the light receptor 
bacteriorhodopsin (bR), organized in trimers and stabi- 
lized by lipids Q. The entire structure appears as a 2D 
hexagonal crystal lattice. The natural role of bR is to 
use sun light for pumping protons outside the cell and in 
doing so, it changes its tertiary structure (conformational 
change). 

Recently, current-voltage (I-V) characteristics of pur- 
ple membrane were analyzed under different experimen- 
tal conditions giving clear evidence of super-Ohmic 
and illumination-dependent responses These mea- 
surements are of paramount importance for the under- 
standing of the mechanism of ET, also in relation with 
the protein activation, furthermore, they are very promis- 
ing for the development of a new generation of organic 
based devices Q • The seminal experiment was carried 
on metal-insulator-metal (MIM) junctions of millimctric 
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diameter, with the insulator constituted by a 5 nm mono- 
layer of PM. The measurements covered a small range of 
bias (0-^1 V), because of the high value of current re- 
sponse (nA level). The response is found to be slightly 
super-Ohmic and grows overall a factor of 2 when the 
sample is irradiated by a green light. These results sug- 
gest that in this protein, like in some organic polymers 
[8|, ET is ruled by tunneling mechanisms. Furthermore, 
the large thickness of PM (5 nm) strongly proposes the 
possibility of multiple carrier jumps across the protein 
(sequential tunneling) S ^• 

In a later experiment [5|, the I-V characterization was 
performed at the nanometric scale, with the technique of 
c-AFM (conductive atomic force microscopy). Accord- 
ingly, one of the contacts is constituted by the tip of the 
c-AFM (100 ~ 200 nm of nominal radius). With respect 
to the first measurement, in the common bias range, the 
current response is lower for about 4 order of magnitude, 
thus the sample is able to sustain higher voltage up to 
about 5 10 V. Measurements were performed without 
any extra light irradiation and with different tip inden- 
tations from about 4.6 nm down to 1.2 nm. At voltages 
above about 2 V, the presence of a cross-over between 
the direct tunneling regime and the injection or Fowler- 
Nordheim tunneling regime was evidenced. 

At present, a microscopic interpretation of the above 
experiments is still in its infancy and, apart from some 
attempts @, Q , a unifying approach able to explain the 
main features and possibly to be predictive for analo- 
gous physical systems is lacking. This paper aims to 
fill this lack of knowledge by implementing a new ET 
model, hereafter called INPA (impedance network pro- 
tein analogue). INPA uses a microscopic description 
of the protein tertiary structure that takes the amino 
acids as single centers of interaction, responsible of charge 
transfer. Finally, the charge transport through a single 
protein results as a process due to the simultaneous ac- 
tivation of multiple pathways in an impedance network 
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on the line of the tunnelling pathway method developed 
by Onuchic and cowkers [lO| . The network approach has 
also strong analogies with a percolative process [11| and 
thus the solution of the transport problem is carried out 
within a stochastic method. With respect to existing 
approaches 12l Il3l|. INPA has some advantages: i) the 



tertiary structure of the protein is directly correlated to 
the macroscopic observables; ii) the different electrical 
responses in the presence/absence of a green light is re- 
produced; Hi) the interpretation of different data [3, S] is 
reconciled iv) it can be applied to other proteins whose 
tertiary structure is known. 

In brief the layout of the INPA is as follows. The pro- 
tein structure is coarse-grained described by means of the 
Cq positions, as obtained by the protein data base (PDB) 
or homology modeling [9*1. Each Cq position is taken as 
corresponding to a node in a graph whose links describe 
the electrical interactions between amino-acids. The de- 
gree of proximity for each node is assigned by a cut-off 
interaction radius, say Rc- In the present analysis, we 
choose Rc = 6 A , a value that optimizes the native to ac- 
tivated state resolution Q. Each link is associated with 
an impedance (a simple resistance in this case) whose 
value depends on the distance between amino-acids. In 
particular: rij = phj/Aij, where p indicates the re- 
sistivity here taken to be same for all the links; p, in 
general, depends on the voltage as detailed below, the 
pedices i,j refer to the amino-acids between which the 
link is stretched, lij is the distance between the labeled 
aminoacids taken as point like centers and Aij is the 
cross-sectional area shared by the labelled amino-acids: 
Ai.j = TT ( — If j/A). The graph is thus mapped into 
an impedance network. To take into account the super- 
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dashed line: Direct tunneling 
continuous line: Direct & FN tunneling 
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FIG. 1: Transmission probability for the reported two tun- 
neling regimes 

linear response, p is chosen to depend on the voltage drop 
as: 

PMAX eV <^ 

piv) = : 

(1) 



where pmax is the resistivity value which should be 
used to fit the I-V characteristic at the lowest voltages, 
Pmin ^ Pmax plays the role of an extremely low se- 
ries resistance, limiting the current at the highest volt- 
ages, and $ is the threshold energy separating the two 
tunneling regimes (a kind of effective height of a tunnel- 
ing barrier). Since ET is here interpreted in terms of a 
sequential tunneling between neighbouring amino-acids, 
the above interpolation formula reflects the different volt- 
age dependence in the prefactor of the current expression 
|15| : I ^ V in the direct tunneling regime, and / ~ 
in the FN tunneling regime. 

For the transmission probability of the tunneling mech- 
anism we take the expression given by Ref . 0, : 
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where Vi,j is the local potential drop between the couple 
of j,j amino-acids and m is the electron effective mass, 
here taken the same of the bare value. 

Figure 1 reports the shape of the tunneling transmis- 
sion probability that includes both the direct and injec- 
tion regimes (continuous curve) together with that cor- 
responding to direct tunneling only (dashed curve). 

In the solution of the resistor network, the tunneling 
mechanism is accounted for by the following procedure. 
First, the network is electrically solved by using the value 
Pmax for all the elemental resistances. Second, by us- 
ing a Monte Carlo acceptance-rejection procedure, each 
Pmax is stochastically replaced by Pmin using the prob- 
ability in Eq. ([2]) and Eq. ([3]) according to the local 
potential drops calculated in the first step. In the high 
voltage region {eVij > <I>), if the stochastic procedure 
gives a rejection, then pmax is replaced by p{Vij) of Eq. 
([T]). The network is then electrically updated with the 
new distributed values of p{Vij). Third, the electrical 
update is iterated (typically 10^ — 10* iterations depend- 
ing on the value of the applied voltage) by repeating the 
second step until the value of the network associated cur- 
rent, converges within an uncertainty less than a few per 
cent. 

To model bacteriorhodopsin in its native state (in 
dark) we have taken the PDB entry 2NTU. The net- 
work is then studied with two different contact configu- 
rations: point- like contacts (in line with previous works), 
and extended contacts. In the first configuration, the net- 
work is connected to the external bias by means of per- 
fectly conductive contacts, the input on the first amino- 
acid, the output on the last amino-acid of the primary 
structure. In the second configuration, the input simu- 
lates the extended tip of the AFM device. Accordingly, 
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FIG. 2: I-V characteristics obtained by simulations with ex- 
tended contacts at the diflerent indentations. The symbols 
and the thin dashed curves refer to calculations. The thick 
continuous line refer to the experimental data in the absence 
of indentation (d=0) when the electrode distance is 4.6 nm ^^l. 
In the inset the theoretical fit for d=0 reported by symbols is 
performed with point-like contacts (see text). 
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FIG. 4: I-V characteristics for the native and activated states 
of bR performed with point-like contacts. Closed and open 
symbols refer to native and activated states, respectively. 
Thick curve to experiments for an electrode distance L=4.6 
nm 0]. In the inset the current is rescaled to reproduce the 
value of the native state at 1 V measured in Ref. (Q)- 
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FIG. 3: I-V characteristics obtained by simulations with ex- 
tended contacts at the different indentations including a leak- 
age current. The thick continuous lines refer to the experi- 
mental data 0] and the symbols to numerical simulations. 

all the nodes with the z-coordinate (direction of the tip 
penetration) larger than that of the first amino acid, as- 
sume the same potential value. The output, point-like, 
remains on the last amino acid. 

The values of pmax , Pmin , and $ are obtained by fit- 
ting the experiments of [5| corresponding to an electrode 
distance of L = 4.6 nm. For this distance, the measure- 
ments are associated with a current crossing a single layer 
of proteins. Actually, we found: pmax = 4 x 10^^ A 
p„j„ = 4 X 10^ n A and 4> = 219 meV. The large differ- 
ence between the pmax and the pmin values is dictated 
by the six order of magnitude spanned by the current 
values. The values of pmax, Pmin correspond to the the 
macroscopic resistance given by the V/I ratio of the ex- 



periment. In such a view, the current measured at 1 V 
in Ref. Q yields a number of trimers equal to about 
N = 10^ for a sample area of 2 x 10^^ nm^, thus to a re- 
sistivity of about 10^° A , for trimer By assuming, 
in present case, the same trimer resistivity, we estimate 
a number of trimers involved in the measured current of 
about N = 10®, thus leading to a crossing area of about 
lO^nm^. However, within a MIM electrical analogue, in 
[H| the effective area deduced by the fit was found to be 
of 0.1 X nm^, about eight order of magnitude smaller 
than what estimated above. This dramatic difference is 
here mainly attributed to the MIM electrical analogue 
used in [Bl that contrasts with the sequential tunneling 
model used here. 

Figure [2] compares the numerical and experimental 
data for the extended contact model and the point-like 
contact model at L=4.6 nm (in the inset). In both the 
cases, the agreement is within the experimental and nu- 
merical uncertainty, and thus considered to be satisfac- 
tory. The threshold energy <i> — 219 meV, is taken to be 
independent of the contact choice. On the other hand, 
when going from the point-like to the extended contact 
configuration the the fitting values of pmax increases 
from 4 X 10^'^ fl A to 8 x 10^^ fl A and also those of pmin 
increases from 4 X 10^ n A to 4 x 10® A. 

The position of the extended contact is then changed 
in order to reproduce the experiments obtained at differ- 
ent indentation of the tip. At increasing depths of the 
extended contact, the net effect is a reduced number of 
amino acids involved in the electrical transport. As a 
consequence, higher currents and a shift to lower poten- 
tial values for the crossover between the direct and the 
FN regime is expected. Numerical calculations reported 
in Fig. [2] confirm these expectations. In particular, the 
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flat contact that simulates the tip indentation was placed 
at the depths of 0.50 nm, and 1.0 nm, from the top of 
the protein. For these depths, the experimental data at 
increasing tip indentations exhibit a quantitative agree- 
ment with the modelling in the region of high voltages 
where FN tunneling prevails. Even if the experimental 
indentation of the tip is greater for about a factor of 3 
with respect to that of the simulations, we consider the 
agreement between theory and experiments to be satis- 
factory in view of the simplifications needed to convert 
the single protein calculation into the macroscopic value 
measured. However, at low voltages the results of the 
simulations underestimate for up to an order of magni- 
tude the values of experiments. The disagreement at low 
voltages is here overcome by assuming the existence of a 
leakage contribution, probably associated with the com- 
plexity of the contact regions, constituted by trimers and 
lipids, with respect to the used model [31 . Accordingly, 
the single resistance associated with each link is replaced 
by the parallel of two resistances, one pertaining to the 
protein and the other to a more realistic modelling of the 
contact regions. 

Figure ^ reports the currents for different indenta- 
tions when an Ohmic leakage current is added to the 
values reported in Fig. 2. The best fit is obtained by 
taking for the leakage resistance the values: 8.4 x 10^^ ft, 
0.11 X 10^2 0.036 X IQi^ respectively for L values of 
4.6, 2.8 and 1.2 nm. The leakage resistance at L=4.6 nm 
is taken equal to the protein resistance value at low bias. 
The decrease of the leakage resistance at increasing in- 
dentation can be related to geometrical effects associated 
with the decreasing of the inter-electrode distance and 
the increasing surface of the tip contact when penetrating 
the protein. Interestingly enough, the microscopic inter- 
pretation is obtained with a length-independent electron 
effective mass, contrary to the case of Ref. [H] where to 
fit experiments the carrier effective mass should increase 
over one order of magnitude at increasing the indenta- 
tion. (We notice that a larger (smaller) value of the 
effective mass would produce a shift to lower (higher) 
value of the transmission probability which can be com- 
pensated by a corresponding change in the value of the 
energy threshold without significant modifications of the 
microscopic interpretation.) Most important, the value 
of the energy threshold is reduced, for about a factor 
of ten, when compared with the values of the barrier 
height found in Ref. Both these features are a con- 
sequence of the sequential tunneling mechanism assumed 
here, and that replaces the single tunneling mechanisms 
of the MIM model previously used ^] . 

As above anticipated, the INPA is able to interpret the 
behaviour of the I-V characteristics carried out with an 
electrode-bilayer- electrode structure and when the pro- 
tein is illuminated or less by green light 0]. To this 
purpose, the model is applied to the PDB entry 2NTW 
(describing the activated state of bR) ^] and the data. 



calculated with the same parameters used for calculating 
the I-V characteristic of the native state, in the extended 
contact configuration. Figure |4] reports the simulated 
data calculated for both the native and activated state. 
The trend evidenced by experiments is here reproduced 
without introducing arbitrary parameters. We remark, 
that present results are compatible with those reported 
by some of the Authors in Ref. Q. Here, the higher value 
of $ used to fit the data of Ref. [1| leads to some minor 
differences in the current responses of the native and ac- 
tivated states which should be justified by the complexity 
of physical system investigated. 

In conclusion, we propose a sequential tunneling mech- 
anism for charge transport in bacteriorhodopsin. The 
model permits a consistent interpretation of a set of ex- 
periments carried out in a wide range of applied electri- 
cal potentials and in the presence or less of an external 
green light. The tertiary structure of the protein enters 
as direct data input and enables one to relate quanti- 
tatively conformational change and sensing action of the 
protein. Finally, also data obtained in very stressing con- 
ditions, like the penetration of an AFM tip in the pro- 
tein membrane, can be finely reproduced. The qualita- 
tive and quantitative agreement between the numerical 
results and experiments poses the INPA, implemented 
here for a sequential tunneling mechanism, as a physical 
plausible model to investigate the electrical properties in 
other proteins pertaining to the transmembrane family. 
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